Nonlinear response of superparamagnets with finite damping: an analytical approach 
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The strongly damping-dependent nonlinear dynamical response of classical superparamagnets is 
investigated by means of an analytical approach. Using rigorous balance equations for the spin 
occupation numbers a simple approximate expression is derived for the nonlinear susceptibility. 
The results are in good agreement with those obtained from the exact (continued-fraction) solution 
of the Fokker-Planck equation. The formula obtained could be of assistance in the modelling of the 
experimental data and the determination of the damping coefficient in superparamagnets. 
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I. INTRODUCTION 

Superparamagnets are nanoscale solids or clusters with 
a large net spin {S ~ 10^-10^). Due to the couphng 
to the environmental degrees of freedom (lattice vibra- 
tions, electromagnetic fields, nuclear spins, conduction 
electrons, etc.), the spin is subjected to thermal fluctu- 
ations and may undergo a Brownian-type rotation, sur- 
mounting the potential barriers created by the magnetic 
anisotropy. This relaxation mechanism was proposed by 
Neel in the late 1940s [H and subsequently reexamined by 
Brown ^] by means of the theory of stochastic processes 
(see also Ref. 0), establishing the basis of the modern 
study of these systems. 

Classical spins with non-axially symmetric Hamiltoni- 
ans can exhibit 3 a large dependence of the thermoac- 
tivation escape rate T on the Landau-Lifshitz damping 
coefficient A in the medium-to-weak damping regime. A 
measures the relative importance of the relaxation and 
the precession in the dynamics. Experiments on individ- 
ual magnetic nanoparticles Q, analysed with accurate 
expressions for the relaxation rate, gave damping coeffi- 
cients in that regime: A ~ 0.05-0.5 [g- 

Uniaxial spins are supposed not to show important 
effects of the damping except in high-frequency condi- 
tions (like FMR experiments). Somewhat unexpectedly, 
non-interacting spins with uniaxial anisotropy, but sub- 
jected to alternate forcing, exhibit a large nonlinear re- 
sponse very sensitive to A 0, which has no analogue in 
the low-frequency linear response. This effect was inter- 
preted in terms of the coupling, via the driving field, of 
the precession of the spin and its thermoactivation over 
the anisotropy barrier. On the other hand, using micro- 
magnetic Langevin simulations 8], Berkov and Gorn 9] 
found that uniaxial spins coupled via dipole-dipolc in- 
teraction also exhibit damping effects such as enhanced 
shifts of the blocking temperature and non-monotonic 
behaviour of the linear susceptibility peaks with the cou- 
pling strength. In Ref. 0| it was shown that these ef- 



fects can be interpreted on the basis of the expression 
derived for T in the mentioned Ref. ;4j , which is valid for 
weak transversal fields but arbitrary damping. Plugging 
heuristically into that relaxation rate the average dipolar 
fields obtained with thermodynamic perturbation theory 
the dynamical effects of the damping on the linear 
response of dipole-dipole coupled systems could be qual- 
itatively reproduced . 

In this article we use a similar analytical approach to 
study the low-frequency nonlinear dynamical response of 
non-interacting classical superparamagnets. We derive 
an approximate expression for the nonlinear susceptibil- 
ity which is in good agreement with the exact (continued- 
fraction) solution of the Fokker-Planck equation. The 
formula obtained is quite simple and may be used to 
model experimental data of the nonlinear response. Ex- 
ploiting its non-trivial damping dependence, the equation 
could assist in obtaining the damping coefficient in these 
systems. The determination of the intrinsic dependences 
of this parameter (on temperature, pressure, etc.) could 
shed some light on the microscopic mechanisms of spin- 
environment coupling in superparamagnets. 



II. BROWN AND KUBO-HASHITSUME 
MODEL 

Let us briefly consider the dynamics of a (sub)system 
accounting for its interaction with the surrounding 
"medium" . This interaction, after the elimination of the 
environmental degrees of freedom, can usually be sepa- 
rated into a time-dependent modulation of the system by 
the proper modes of the environment (fluctuating term), 
and the back-reaction on the system of its action on the 
surrounding medium (relaxation or damping term). 

This approach was particularized phenomenologically 
by Brown 2J and Kubo and Hashitsumc 3] to classi- 
cal spins, by introducing the stochastic partner of the 
Landau-Lifshitz equation. The associated Fokker-Planck 
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equation governing the time evolution of the probabiHty 
density of spin orientations P{s) can be written as |l2j| 

2TDdtP^iJ-[{B,-{sAB,)+ij\P . (1) 

Here = —is A {d/ds) is the generator of infinitesi- 
mal rotations and = —f3{d'H/ds) is an effective field. 
The Landau-Lifshitz relaxation parameter A (dimension- 
less) measures the relative importance of the damping 
and precession terms. Finally, td is the relaxation time 
of isotropic spins (the counterpart of the Debye time in 
dielectrics) 



1 

TD = T 



A 2jkBT ' 



(2) 



where m is the spin magnitude and 7 the gyromag- 
netic ratio. For generalisations of the Brown and Kubo- 
Hashitsume model, see, for instance, Refs. 



Equating the coefficients with the same oscillating factor 
exp{kiu!t) one gets the ATV^, which are directly related 
with the susceptibilities. We keep at each order k only 
the leading term in (if required, the next order terms 
can be obtained along the same lines). 

Let us assume that in the absence of perturbation the 
two states (wells) are equivalent (symmetric). Taking 
into account that {A~^ + A^) is a relaxation rate (and 
hence even in ^) and that (^+ — A~) is related to the 
forcing (odd in ^), the response will depend only on the 
sum of the A^ +A^ for even k and the difference A'^—A'^ 
for odd k [vd. Eq. Taking this into account, the 

amplitudes of the response read (AA'o = and AiV2 — 0) 



AiVi = - 



Fq + 'ilUJ 



(Fo + lu) (Fo + 3ic^) 



III. GENERIC BALANCE EQUATIONS 

Before proceeding from the Fokker-Planck equation 
to study the nonlinear dynamics, let us consider some 
generic expressions for systems describable by a set of 
kinetic balance equations for some occupation numbers 
7V+ and 7V_ 



where we have introduced the relaxation rate in the ab- 
sence of forcing 



Fo = + A^ 



(9) 



These results are quite generic. In particular cases the 
A^ will be constructed from the especific details of the 
model. 



N+ = -A+N+ + A-N^ 
iV_ = A+N, - A-N^ 



(3) 



IV. BALANCE EQUATIONS: SPIN DYNAMICS 



Here the A^ are some transition amplitudes which de- 
pend on the external forcing or control parameter ^. The 
occupation numbers satisfy the "constrain" — 1, 

which indicates the conservation of the number of rep- 
resentative points (iV+ = — iV_). The response of the 
system is characterised by the difference in populations 
A7V = N+ - N^, which obeys 



dt 



-AN 



(A^ 



A-] AN - A+ - A- 



(4) 



Thus, {A~^ +A~) plays the role of a relaxation rate while 
the inhomogeneous term (A+ — A~) is to be related to 
the external forcing. 

To get the linear and first nonlinear susceptibilities 
(or corrections to the linear susceptibility due to a weak 
static forcing) we expand A^ in a series of powers of ^ to 
the third order 



A' 



^0 



^Af+eAt+eAi 



(5) 



Let us consider in detail the case of harmonic forcing 
A^{t) = iA^(e+''^*-he-'"*). First, we replace ^ by A£,{t) 
in the above expansion. Next, we plug into the dynamical 
equation for AiV both the resulting A^ (t) and the Fourier 
expansion of the population difference, namely. 



AN ~ AiVo + (^)AiVie" 



(^)'A7V2e2-* 



\ 2 J 



{^fAN.e'^'^ 



(6) 



For a spin with the simplest uniaxial anisotropy in a 
field (chosen by convenience to lay in the XZ plane), the 
Hamiltonian can be written as (s = rn/m) 



'PH = asl+£,«s,+^A 



(10) 



The anisotropy term has two minima at S2 = ±1 (the 
"poles") with a barrier between them at Sz = (the 
"equator"). The spin-Hamiltonian parameters are in- 
troduced in temperature units: a =^ D/k^T is the 
anisotropy barrier while and are the longitudinal 
and transverse components of the field ^ = mB/k^T, 
with respect to the anisotropy axis. 



A. Balance equations 

Garanin et al. 4] rigorously derived from the Fokker- 
Planck equation a set of balance equations for the occu- 
pation numbers in the upper Sz > well (our 7V+) and 
the lower Sz < well (-/V_), namely 



F (iV;^iV_ - Nl''N+) 



7V_ = -F {Nl'^N^ - Nl'^N+] 



(11) 



Here N^ = 2±/Z are the equilibrium occupation num- 
bers with Z± the partition function restricted to the up- 
per and lower wells, respectively. On comparing with 
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the generic Eq. (jj)), we find for the transition amphtudes The relaxation rate F is given by 
(note the sign reversal) 

A± = TN^"^ . (12) 



1 



kBT 
m/7 



2tt 



dip e 



ds. 



(13) 



where Sz and ip (the azimuthal angle) are the canonical variables of the spin. The function Ci^z, v) is determined by 
the quasistationary solution of the Fokker-Planck equation 



dip dsz dsz dp 



dn , d 
h kBT— 

OSz OSz 



1 



dC 



1 



1-si 



dp dp J dp> 



subjected to the boundary conditions C = 1 s-^d C — i'^ th e bottom of the lower and upper wells, respectively. 



B. The transition amplitudes at low fields 

In order to get the field expansion of the transition 
amplitudes = TN^, we need the corresponding ex- 
pansions of the equilibrium occupation numbers and the 
relaxation rate. 



1. The low-field equilibrium occupation numbers 

The partition function corresponding to the Hamilto- 
nian IjlUfl can be written as 

Z = y^ds,^ ^exp{~m) • (14) 

The one-well partition functions Z± correspond to in- 
tegrate Sz over [0,1] and [—1,0], respectively. Writting 
Sx = y^l — Sz COS ip and following Shcherbakova '15] in 
doing first the integrals over p, we get the unified ex- 
pression 

Z±= [ dsz exp(as2 ±^||s,)/o(e±Vl - , (15) 

where we have made the change of variable Sz — > —Sz 
in Z_ and Iq is the modified Bessel function of the first 
kind of order lOj Sec. 11.5]. Thus, Z± can be written in 
terms of an integral over Sz only [cf. Eq. (11) in Ref. p^]. 

Calling a = i^\\Sz and b = ^j_-\/l — and using the 
expansion (to third order) Io{b) ~ 1 + ^b^ , we obtain for 
the field-dependent part of the integrand of Z± 

e"/o(6) ~ 1 + a + ia^ + + ia^ + ^^b^ . (16) 

Now we introduce the zero-field averages in one well 

(4)w = Y,""; , , (17) 



where the denominator is proportional to the zero-field 
partition function Zq = J -^dsz e'^*^ . Introducing the ex- 
pansion of e°/o(^) in Eq. H15|l we get the field expansions 
of Z± and, by adding them, that of Z itself. Using the 
binomial formula to get the corresponding expansion of 
1/Z, and multiplying this by those of Z±, we finally ob- 
tain the equilibrium occupation numbers = Z±/Z. 
These can be written as (note that N'^ + N'^ = 1) 

7V^'>=i(l±Zie||± 23^11 ±Zl,2e||Ci) , (18) 

with the coefficients Zj^^ given by 

zi = (s,)w (19) 
Z3 = H(^'>w-3(s.)w(s^)w) (20) 
^1,2 = -3((^')w-(s.)w(s^)w) . (21) 

In the coefficients Zj.^ the first index is the power of ^|| 
and the second (omitted when zero) the power of 
Note that the expressions for are valid for an arbi- 
trary uniaxial potential. 

2. The low-field relaxation rate 

As the spins have inversion symmetry in the absence 
of the field, the total relaxation rate should be an even 
function of the field (F accounts for jumps over the en- 
ergy barrier in both directions). For spins with uniaxial 
anisotropy we can write 

F~Fo(l-K5||^f + <?±^i) , (22) 

where Fq is the zero-field relaxation rate and the expan- 
sion is valid to third order. The vanishing of the term 
c< follows from the invariance of the relaxation rate 
upon field refiection through the barrier plane in uniaxial 
spins. 
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3. Generic expression for the transition amplitudes 

Plugging the expansions ((T5)l and in the expression 
for the relaxation amplitudes ~ TN'^, we arrive at 



± 



(^3 + zig\\) + (zi,2 + zig±) 6||6 



(23) 



where we have introduced the direction cosines of the 
field 



Let us write the components A'^ ± A^ that enter the 
equations for the response [Eqs. Q"®]- Note first that 
the obtained result fulfils the consequences of the well- 
symmetry mentioned above A'^ — A^ = and A2 — A2 — 
0, along with A'^ + A^ =0 and A^ + A^ = 0. The 
combinations entering in the response are given by 



A^r, 



-(i+-ir) = ^1^11 



At 



A-, 



9ib\ 



A-, 



(z3 + zi^ii) h\ + (zi,2 + ^i.g±) 



where we have introduced the notation A — A/Fo. 



V. EXPRESSIONS FOR THE DYNAMICAL 
SUSCEPTIBILITIES 

A. Generic expressions 

Let us divide numerator and denominator of the AA^^ 
by To and introduce the relaxation time 



(24) 



Then, Eqs. ((7J) and ||SJ| appear with ii's (= AjV^) in the 
numerator and factors 1 + fciwr in the denominator. The 
response is the projection of the average spin onto the 
field direction. This projection is obtained by multiplying 
the difference in the populations of both wells AA^ by 6|| 
and the spin magnitude m. 

In order to get the susceptibilities, recall that we used 
the field in temperature units, ^ = mB/k^T, which 
yields factors (m/fceT)*^. Thus, x^^'^ = {m/k^Tf/^Nk x 
(m6||), and the linear susceptibility reads 



X 



(1) 



m b 



2 f,2 

I 



zi 



knT 1 



(25) 



For the first nonlinear susceptibility we obtain 

(3) _ "^^ ^1 Z3 + (Z3 + ziffii) iujT 



(fcBT)3 (1 + iut) (1 + Silut) 

{kBTY (1 + iwr) (1 + Siwr) 



(26) 



where we have grouped terms with the same powers of 
so that the angular dependence (tensor structure) 
is better recognised. 

Note that the expressions for the response are quite 
generic and depend only on the coefficients of the ex- 
pansion of the equilibrium occupation numbers and the 
relaxation rate. Specific formulae will be obtained de- 
pending on the features of the uniaxial potential and the 
approximations done in calculating the coefficients Zj^^ 
and 



B. The case of low temperatures 

Let us now specialise the above formulae to the case of 
low temperatures, where the superparamagnetic blocking 
takes place for long measurement times (or equivalently 
low frequencies, as those of ordinary magnetic experi- 
ments) . 

The coefficients Zj ^^ are determined by the one-well 
averages of low-order powers of Sz [Eqs. H19(l - (|21|l ]. For 
anisotropy energy cx , these can be obtained along the 
lines of the calculation of Ref. [11 App. A]. Thus, us- 
ing the asymptotic expansion of the confluent hypergeo- 
metric (Kummer) functions [161] . one finds the following 
low-T expansion 



4ct2 



(27) 



where a = D/k^T. With this result, we immediately 
obtain the coefficients of the field-expansion of the equi- 
librium occupation numbers 



zi 
z-i 
z\a 



1 



1 

2^ 



1 

3 
1 

8^ 



1 

2^ 



1 

2^ 

1 

^ 4^ 



(28) 
(29) 
(30) 



For the relaxation rates we shall only obtain the leading 
order term in 1/ct; for consistency, the above Zj^^ will 
only be used up to such order (we shall return to this 
point below). 

To get the coefficients g\\ and appearing on the field- 
expansion of the relaxation rate [Eq. H22|l ] , one can choose 
special configurations in which they are known (strictly 
longitudinal and transversal fields) |lfl| |. Expanding the 
formula for F in the presence of a longitudinal field ^|| 
Ulli one finds 



F(e||,e± = 0)^Fo(l + ieff) , Fo 



1 2 



TD 



3/2 - 

a ' e 



where Fq is Brown's zero-field result for the relaxation 
rate. Comparison with the general expansion H22|l gives 
the longitudinal coefficient = 1/2. Note that in a 
longitudinal-field the damping parameter A only enters 
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through td [Eq. (01] and hence A only matters to es- 
tabUsh a global time scale. In other words, results for 
different damping parameters presented in units of td 
show complete dynamical scaling, and in this sense the 
A dependence is said to be trivial. 

Non-trivial effects of the damping arise in a transversal 
field. Nevertheless, there is no general expression for the 
relaxation time in the presence of transverse fields. In 
Ref. however, a low-temperature formula valid for 
weak transversal fields was derived, which is perfectly 
suited for the purpose of determining g^, namely 



r(6| =0,6 



1 



= Aai/2. (31) 



Comparing with the expansion (|22(l one gets the transver- 
sal coefficient g± = F/4. 

Summarising, the low-temperature coefRcients in the 
field expansion of the relaxation rate [Eq. H22(l ] are given 
by 



9\\ 



9± 



(33) 



The function F{a) takes into account, without further 
approximations, the effects of the damping and is given 
in terms of the incomplete gamma function 7(0; z) = 
dtr^^e-* by 



F{a) = 1 + 2(2a^e) 



2pU/(2a') 



7 1 + 



1 



1 



2a2 ' 2a2 



(32) 



The function F decreases towards 1 for strong damping 
[vd. Eq. (|1T|| below]. Then, and g±^ are of the same 
order of magnitude. However, F grows as 1/A for weak 
damping, where the relaxation time turns to be very sen- 
sitive to the damping (or to transversal fields). 



Before giving the corresponding expressions for the suceptibilities, we write the complete formulae for the transition 
amplitudes — TN'^ at low temperatures (including only the leading order in the l/cr expansions) 



1±6||6 



(34) 



This expression can also be of use in other problems like the obtaining of the linear susceptibility in weak bias fields. 
The linear susceptibility arising from Eq. (|25|l simply reads 



X 



(1) 



k^T 1 -f lOJT 

while the nonlinear susceptibility obtained from Eq. H26|) is given by 



X 



(3)^_1 - 



1 



3(fcBr)3 (l + icjT)(l + 3iwT) 4 (fcBT)3 



F(a) 



(1 + iwr) (1 + Siwr) 



(35) 



(36) 



This is one of the main results of this article. As special 
cases we consider the response to a strict longitudinal 
field (&11 = 1 and = 0) 



This particular case (with I/ct corrections) was derived 
in Ref. 20] from an analytical treatment of the Fokker- 
Planck equation disregarding the precession terms. 



X 



(3) 



(37) 



3 [kBTf (1 + it^T) (1 -|-3iwr) ' 
and the response of an ensemble of identical spins with 



axes distributed at random (fejj 
1 1- 



X 



(3) 



: 1/5 and = 2/15) 
(38) 



15 (fcBT)3 (1 -I- ilot) (1 -f Siwr) 



Recall now that for overdamped spins one has F ^ I. 
Therefore, the formula derived reduces in this case to 



^ Ia»i 



1 



1 



15 (fcBr)3 (I + ilut) (1 -|-3iwr) 



(39) 



C. Practicalities implementing the analytical 
expressions 

To compensate for disregarding l/cr corrections, we 
can heuristically replace the low-T equilibrium linear 
and nonlinear susceptibilities by their exact expressions 
18, 21, 22]. Further, for axes at random, we can sim- 
ply use the leading correction for the equilibrium non- 
susceptibility 



(3) 
Xoq 



1 



15 {kBT} 



(40) 
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As for the function F{a), we can use the approximate 
forms 4] 



F{a) 



7r/a — 1/3 + \/7ra/6 , a<l 



(41) 



FinaUy, for the relaxation time at zero field, we use the 
accurate in terp olation formula of Gregg, Crothers and 
Wickstead IM 



To 



1 



t3/2 



td V 1 + l/o- 



a2- 



1 



exp((T) — 1 



(42) 



With these prescriptions, no special functions appear in 
the analytical formulae, which are expressed in terms of 
simple functions and polynomials. We shall see below 
that the agreement of the resulting formulae with the 
exact continued-fraction results is quite good. 



VI. RESULTS FOR THE NONLINEAR 
SUSCEPTIBILITY 

Numerically exact results for the nonlinear susceptibil- 
ity can be obtained by solving the Fokker-Planck equa- 
tion by continued-fraction methods. In this approach, 
one considers the equations for the spherical harmonics 
Y^^Szt'^) averaged with respect to the non-equilibrium 
distribution Pis) obeying Eq. QJ. The equations for 
the 17" 0, 125I (see Ref. [I^i for an alternative deriva- 
tion) can be solved perturbatively in the forcing A^(t) 0- 
At each perturbative leveL on introducing apropriate 2- 
vectors and 2x2-matrices 0,123 the equations for the 
can be cast into the form of a three-term recurrence rela- 
tion (in the index i with fixed m). This recurrence can be 
solved efficiently and accurately by using matrix contin- 
ued fraction methods |2^. Finally, the average response 
of the system is obtained with help from the relations 
Sz = Y^ and + isy = -^8^3 ^1 ■ 

The features of the nonlinear susceptibility spectra of 
classical superparamagnets in the experimentally most 
common case of anisotropy axes distributed at random 
are displayed in Fig. ^ The nonlinear susceptibility x''^'' 
shows a large dependence on A 0, dependence that is 
absent in the linear susceptibility for the same axes dis- 
tribution and in the longitudinal and strict transverse 
nonlinear susceptibilities. The sensitivity to the damping 
was interpreted in terms of the dynamical saddle point 
created by the oblique driving field in the uniaxial poten- 
tial barrier. This saddle favors inter-potential- well jumps 
that would be unlikely if the field were in the linear range 
(weakly deformed barrier), and hence leads to an increase 
of the magnitude of the low-T response. 

To illustrate this interpretation, consider one spin that 
after a "favourable" sequence of fluctuations, approaches 
the top of the barrier but does not surmount it. In the 
subsequent spiralling down back to the bottom of the 
well, a strongly damped spin descends almost straightly. 
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FIG. 1: Nonlinear susceptibility x'^' of classical spins with 
randomly distributed anisotropy axes vs. frequency. The tem- 
perature is fcer/D = 0.05 (cr = 20) and A = 1, 0.1, 0.03, and 
0.01 (bottom to top). The lines are the approximate formula 

rand the symbols the exact continued-fraction solution 
The upper and lower panels display, respectively, the real 
and imaginary parts of x'^'('^). 



whereas a weakly damped spin executes several rota- 
tions 1/'^) about the anisotropy axis. This allows 
the weakly damped spin to pass close to the saddle area, 
where it will have additional opportunities, not available 
for the damped spin, to cross the barrier, enhancing its 
relaxation rate. Naturally, this mechanism will make a 
difference at low temperatures, where reaching the bar- 
rier region is a rare event. 

The analytical expression derived in Sec.^is displayed 
for comparion in Fig. ^ The agreement with the exact 
results is notable. The maximum relative error occurs 
when A = 0.01 at the peak of x*-'^-''(^^), and it is only a 
2%; going to much higher temperatures, a — 10, where 
the approach should start to fail, that error is still less 
than a 4% (recall that typical experimental conditions 
correspond to tr ~ 20-25). This agreement, in turn, sup- 
ports the interpretation discussed above of the damping 
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dependence. The reason is that our analytical expres- 
sion includes at its heart the formula for F in a weak 
transverse field Q, which accounts for the effects of the 
corresponding saddle point on the relaxation rate. 

VII. DISCUSSION 

With help from the expression for the relaxation rate 
of Garanin et al. jj], it has been shown that effects ap- 
parently different as the damping dependence of F in 
transverse fields, the damping sensitivity of the nonlinear 
susceptibility , and the effects of the damping on the 
linear response of dipole-dipole coupled systems 0, ^3 j 
have a common origin, which is the strong sensitivity of 
the relaxation rate to the damping in transverse fields. 

The damping parameter A carries information of the 
coupling with the environmental degrees of freedom caus- 
ing the relaxation in superparamagenets. In Ref. it 
was suggested that the large A-dependence of x*^^^ could 
be exploited to determine A experimentally. This would 
by-pass its determination from the pre-exponential factor 
To in the relaxation time (cx td), which usually only gives 
an order-of-magnitude estimate. Besides, one does not 
need high frequencies to explore the effects of A, in con- 
trast to magnetic resonance experiments (avoiding the 
associated technical difficulties). Clearly, the availability 
of a simple analytical expression to model the x^'^' data 
should be of great assistance to determine A. 

Coffey et al. j22ij suggested that a method based on the 



linear susceptibility with superimposed bias fields would 
made the resort to the nonlinear response unnecessary. 
However, in their case the fittings should be done to a 
formula involving more complicated expressions, both for 
the equilibrium parts (numerically obtained) and the re- 
laxation times (involving the different Kramers' regimes). 
In addition, their results have to be eventually integrated 
over the distribution of anisotropy axis orientations. 

Our expression is free from those complications (basi- 
cally involves the simple zero-field r of Brown and equi- 
librium susceptitibilities known analytically). Besides, 
the measurement of x^'^-' is becoming standard nowadays 
(see, for instance, Ref. "2^). In the modelling of real 
experiments the incorporation of the particle-size distri- 
bution can be done by simple integration of our equa- 
tion. For these reasons, we consider the method based 
on x*''^^ more suited for the experimental determination 
of the important, and hitherto quite evasive, dissipation 
parameter in superparamagnets. Finally, the genericity 
of the intermediate expressions derived could allow the 
incorporation of quantum effects, by taking them into 
account in the field-expansion coefficients of the equilib- 
rium quantities and the relaxation rate. 

Acknowledgments 

This work was partially supported by DGES (Spain), 
project BFM2002-00113 and the Swedish Foundation for 
Strategic Research (SSF). Discussions with P. Svedlindh, 
P. Jonsson, and F. Luis are warmly acknowledged. 



[1] L. Neel, Ann. Geophys. 5, 99 (1949). 
[2] W. F. Brown, Jr., Phys. Rev. 130, 1677 (1963). [16 
[3] R. Kubo and N. Hashitsume, Prog. Theor. Phys. Suppl. 

46, 210 (1970). [17 
[4] D. A. Garanin, E. C. Kennedy, D. S. F. Crothers, and 

W. T. Coffey, Phys. Rev. E 60, 6499 (1999). [18 
[5] W. Wernsdorfer et al, Phys. Rev. Lett. 78, 1791 (1997). [19 
[6] W. T. Coffey et al, Phys. Rev. Lett. 80, 5655 (1998). [2Q 
[7] J. L. Garcia-Palacios and P. Svedlindh, Phys. Rev. Lett. 

85, 3724 (2000). [21 
[8] A. Lyberatos, D. V. Berkov, and R. W. ChantreU, 

J. Phys.: Cond. Matter 5, 8911 (1993); J. L. Garcia- [22; 
Palacios and F. J. Lazaro, Phys. Rev. B 58, 14 937 
(1998). [23 
[9] D. V. Berkov and N. L. Corn, J. Phys.: Condens. Matter 

13, 9369 (2001). [24 
[10] P. E. Jonsson and J. L. Garci'a-Palacios, Europhys. Lett. 

55, 418 (2001). [25 
[11] P. E. Jonsson and J. L. Garcia-Palacios, Phys. Rev. B 

64, 174416 (11 pages) (2001). [26 
[12] J. L. Garcia-Palacios and P. Svedhndh, Phys. Rev. B 63, 

172417 (4 pages) (2001). [27 
[13] D. A. Garanin, V. V. Ishchenko, and L. V. Panina, Teor. 

Mat. Fiz. 82, 242 (1990), [Theor. Math. Phys. (USSR) [28 
82, 169 (1990)]. 
[14] J. L. Garci'a-Palacios, Eur. Phys. J. B 11, 293 (1999). 
[15] V. V. Shcherbakova, Izvestiya, Earth Phys. 14, 308 



(1978). 

G. Arfken, Mathematical Methods for Physicists (Aca- 
demic Press, Boston, 1985), 3rd ed. 

P. J. Gregg and L. Bessais, J. Magn. Magn. Mater. 203, 
265 (1999). 

J. L. Garci'a-Palacios, Adv. Chem. Phys. 112, 1 (2000). 
A. Aharoni, Phys. Rev. 177, 793 (1969). 
Y. L. Raikher and V. L Stefanov, Phys. Rev. B 66, 
214406 (2002). 

J. L. Garci'a-Palacios and F. J. Lazaro, Phys. Rev. B 55, 
1006 (1997). 

J. L. Garci'a-Palacios, P. Jonsson, and P. Svedhndh, 
Phys. Rev. B 61, 6726 (2000). 

P. J. Gregg, D. S. F. Crothers, and A. W. Wickstead, 
J. Appl. Phys. 76, 4900 (1994). 

Y. P. Kalmykov and W. T. Coffey, Phys. Rev. B 56, 3325 
(1997). 

Y. P. Kalmykov and S. V. Titov, Phys. Rev. Lett. 82, 
2967 (1999)'. 

H. Risken, The Fokker-Planck Equation (Springer, 
Berlin, 1989), 2nd ed. 

W. T. Coffey, D. S. F. Crothers, Y. P. Kalmykov, and 
S. V. Titov, Phys. Rev. B 64, 012441 (4 pages) (2001). 
P. Jonsson, T. Jonsson, J. L. Garcia-Palacios, and 
P. Svedlindh, J. Magn. Magn. Mater. 222, 219 (2000). 



